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. Abstract 

We perform numerical experiments to study the Lyapunov spectra of dynamical 
systems associated with the Navier-Stokes (NS) equation in two spatial dimensions 
truncated over the Fourier basis. Recently new equations, called GNS equations, 
have been introduced and conjectured to be equivalent to the NS equations at large 
ff^ . Reynolds numbers. The Lyapunov spectra of the NS and of the corresponding GNS 

\ systems overlap, adding evidence in favor of the conjectured equivalence already 

^ I studied and partially extended in previous papers. We make use of the Lyapunov 

CN ■ spectra to study a fluctuation relation which had been proposed to extend the 

, "fluctuation theorem" to strongly dissipative systems. Preliminary results towards 

the formulation of a local version of the fluctuation formula are also presented. 
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1 Introduction 



The idea of viewing turbulent fluid motions as trajectories in the phase space 
of a chaotic dynamical system is a long standing one. In particular, treating 
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turbulence in terms of invariant probability distributions on the phase space 
of the dynamical system given by the Navier-Stokes (NS) equation for a fluid, 
and asserting that the distributions ought to have special properties, is an 
idea which dates back to the early 1960 's: it has been formalized a little later, 
see[l].Q 

Mathematical and numerical studies have been performed, at first for the 
purpose of relating fluid dynamical parameters, such as the Reynolds number 
or the amplitude of external forcing, to the fractal dimension of the attractorQ 
In this way the onset of turbulence has been quite well understood, along with 
the discovery that the relevant physical phenomena are far more interesting 
and complex than previously anticipated, see[3,4,5,6,7]. 

However, understanding developed turbulence remains a most challenging 
problem and it is not clear how to use the basic ideas of Lorenz, Ruelle and 
Takens in this context. The approach to the analysis is very often based on 
numerical experiments in which simulations of reduced models of turbulence 
are considered. For instance the shell models of Yamada and Okhitani [8,9,10] 
have been studied in great detail with the aim of finding or checking connec- 
tions between Lyapunov spectra and the turbulent cascade. 

Another viewpoint [11,12] stems from the well known fact that macroscopic 
irreversibility of particle systems can be obtained from microscopic reversible 
dynamics, provided the ratios between the microscopic and the macroscopic 

^ Quoting [2], "It is widely accepted that the statistical properties of a turbulent 
flow can be described by a dynamical evolution on a strange attractor in a high- 
dimensional phase space. An important characterization of the ergodic properties 
of the attractor is given by the Lyapunov spectrum." 

^ It is important, when discussing fractality, to distinguish between attractor and 
attracting set. An attracting set is a closed set such that all points in its vicinity 
get as close as one wants to the set as time tends to infinity and, furthermore, such 
that none of its closed subsets have the same property. The notion of attractor, 
that none of its closed subsets have the same property. The notion of attractor is 
associated with a method of choice of initial data in the vicinity of an attracting 
set. Suppose that the data are randomly selected with a distribution /iq and, with 
/io^probability 1, their asymptotic motion is described by a probability distribution 
fi on the attracting set, in the sense that time averages can be obtained simply 
by integration with respect to //. Then, any subset of the attracting set with fi- 
probability 1 and smallest Hausdorff dimension is called "attractor" for the random 
choices with distribution fiQ, and is called the statistics of fiQ. Quite commonly, the 
attracting set is not fractal, and it might even coincide with the whole phase space, 
while the attractor for (some) //q is fractal and with dimension lower than that of 
the attracting set. We recall that the special statistics generated on an attracting set 
from data randomly chosen with a probability distribution /iQ absolutely continuous 
with respect to the volume is called a "iSiJ-B distribution" , which is usually unique: 
in this paper we shall, as usual, only consider the SRB statistics. 
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time and length scales are sufficiently large. In that case, in fact, the micro- 
scopic reversible equations of motion are "equivalent" to dissipative equations 
of motion at the macroscopic level. 

Assuming that forced turbulence is described by the (irreversible) NS equa- 
tion, a different (reversible) equation called the "Gaussian Navier Stokes", 
or GNS, equation was proposed to describe the same physical phenomena 
[11,12]. In the GNS equation the usual viscosity constant is replaced by a 
fluctuating term, which makes the equation itself time reversal invariant and 
keeps the total energy exactly constant. Equivalence should occur when the 
constant energy of the GNS evolution equals the average energy in the NS 
evolution provided the Reynolds number is large: in a way which is similar to 
the equivalence between microcanonical and canonical ensembles in statisti- 
cal mechanics [13], with the large volume limit replaced by the limit of large 
Reynolds number R. The conjecture states that the forward time statistics 
of the NS and of the GNS models of the same fluid are the same, for what 
concerns "local observables" , i.e. for the time averages of functions of the ve- 
locity field which depend only on a small number of its Fourier components 
physically corresponding to observables on the inertial scale. 

The idea, if correct, could be used to analyze phenomena which would not 
be easily understood in terms of the NS equation. For instance, time reversal 
invariance and chaoticity of the GNS dynamics imply that the fluctuations 
of certain quantities verify a relation (cf. Ref.[14]), here referred to as "fluc- 
tuation relation" (FR), which might eventually be amenable to experimental 
verification as attempted in Refs. [15,16]. It should be remarked that the pos- 
sibility of using several different representations of the statistical properties of 
a turbulent fluid flow, which result equivalent from a practical point of view, 
was considered already in previous works, see, e.g. Refs. [17, 18, 19], and is still 
under investigation [20]. 

Of course despite the expected equivalence, the NS and GNS representations 
of a turbulent fluid remain different, similarly to the microcanonical and the 
canonical ensembles descriptions of equilibrium statistical mechanics: although 
"equivalent in the thermodynamic limit" they differ in many respects {e.g. in 
the distribution of fluctuations of the total energy) . 



2 NS and GNS equations. Aims of our analysis. 

Suppose that the side of the periodic cell is L, and that the forcing term f 
is the product of a dimensional parameter F times a flxed function. We shall 
suppose the Fourier components Fk to have modulus 1 and, usually, to vanish 
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for all but one value of k.[^ Then, in the NS case with viscosity u, we define 
the Reynolds number R to be 

= FL^iy-^, or R = ly-^VWI^ . (1) 



and we shall consider the NS equations in dimensionless form. Here, the lengths 
are measured in units of L/27r and time in units of L'^/iy, so that the change 
of variables u{x,t) = yu9(x/L,Ct) with V = FL'^/p,C = i^/L"^, yields (re- 
defining simply by u the dimensionless field for notational convenience) the 
dimensionless NS equations 

u + R^{u-d)u = Au + i -dp , d-u = (2) 



with u periodic in the position coordinates with period 27r and divergenceless 
{i.e. u represents an incompressible flow). Correspondingly, the GNS equations 
will be deflned by 

u + R^{u-d)u = aAu + i-dp , d-u = Q (3) 



where i? is a free parameter and a is deflned so that the solutions of the 
equation (3) conserve Qo = J u^rfx, i.e. twice the total kinetic energy. 

In order to compare motions that the GNS equations generate with (suitably 
selected) motions of the NS fluid, the constant Qq has to be chosen. For a 
given R and taking as reference the NS equation, we call (Qo) the inflnite time 
average of Qq produced by the NS evolution. We shall compare, according to 
the analysis of [21,22], the solutions of the NS equation with those of the GNS 
equations with the conserved quantity Qo set equal to {Qo). 

We write u(x) = J2k e*'^'^Wkk-'-/|k|, where the scalar term obeys = —uZ^, 
and k"*- = {k2, — fci) if k = {ki, k2). Hence Qo = (2vr)2 X^k [""kP (and the total 
vorticity is Qi = {^nYJ^k |kp|MkP)- If we deflne the ejfort Gi as 



Gi 



(2vr)^E 



a — f — dp) , (a — f — dp) 

ilk + iR^ Ej+h=k ^^^^ijhpp^j^h - /k)' 

k2 



(4) 



where (■, ■) denotes the usual scalar product in L2{dx.), and a = du/dt = 
dtu + R'^{u- d)u is the acceleration fleld, Eq.(3) can be obtained from Gauss' 
least effort principle by considering an Euler fluid subject to the constraint 

^ The cases of Fig. 2 are the only exceptions considered here. 
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of constant Qq, i.e. constant total energy. In this case the acceleration a of 
the velocity field u minimizes the effort Gi among all possible accelerations 
compatible with the incompressibility and with the constraint. If we write 
Eq. (3) in terms of Fourier components we find 

u^^-iR" Yl nil. 11.1 ^ ^j^h-«k^t^k + /k (5) 
j+h=k 



where a is given by 



Ek /kl^k /„x 



Of course the fact that the equations above can be obtained from a variational 
principle might be merely accidental, and in itself does not justify their use. 
However in Refs.[ll,21,22,12] it has been argued that they should give a de- 
scription of the fluid motions which in the large Reynolds number regime is as 
accurate as the NS equations themselves. More generally, one can consider an 
Euler fluid subject to a Gaussian constraint on with effort function G^, 
£ < m + 1, given by 

g^ = (27r)2^|k|2-|i.kr, 

k 

G,,^ = ((A)-(^-'")(a-f-ap),(a-f-9p)) . (7) 

This leads to equations like (3) with the Laplace operator replaced by A^, 
and with a replaced by a new coefficient Pi^^ thai can be easily computed by 
imposing that Qm remains exactly constant. We call the resulting equations 
" G-hyperviscou^' equations while we call "hyperviscous" the corresponding 
equations in which is replaced by a constant. The above case (5) corre- 
sponds to £ = 1, m = 0. The special cases i = 0, m = and £ = 1, m = 1 have 
been considered in [11] and the general hierarchy of the hypcrviscous equa- 
tions has been considered and partially studied in [22]. The case i = l,m = 1, 
for instance, yields: 



Ek 


|k| 


P/k^k 


Ek 


|k| 


4 


Mk 


2 



= ^^.l: . ■ (8) 



A general feature of the G-hyperviscous equations is that they generate re- 
versible dynamics in the sense that the velocity reversal operation lu^ — — tik 
anticommutes with the time evolution map St (the mapping of an initial da- 
tum u into its value S'^u at time t) in the sense that 

I-St^S.fl . (9) 



5 



Note that here we call GNS equations the case i = l,m = for consistency 
with [22] because we shall often refer to Ref.[22], although in [11] this name 
was used for the case £ = l,m = 1. 

In this paper we study via computer simulations: 

• (1) Relations between the Lyapunov spectra of corresponding NS and GNS 
models. Equivalence between the Lyapunov spectra is not part of the earlier 
proposals on equivalence because the Lyapunov exponents are not "local 
observables" : however it seems possible that they reflect properties of the 
individual modes of the velocity field [23] so that it is interesting to test 
whether equivalence extends to them. In the models that we consider the 
Lyapunov spectra appear to be essentially identical within the errors. 

• (2) Fluctuation properties of the observable defined by the phase space 
contraction in the reversible GNS equation and its relation with the Lya- 
punov spectra following the proposal in [11]. The relation proposed in [11] 
is, however, based on several assumptions which we cannot check directly. 
Therefore we have measured quantities for which predictions are available 
as consequences of the mentioned assumptions, finding that the predictions 
match with the observations. This lends support to the idea that the prop- 
erty called axiom C" as well as the assumptions on the structure of the 
attracting sets discussed in [11,24] may be actually verified in the models 
considered here. 

• (3) The equivalence analysis can be extended to involve not only NS and 
GNS equations but also to several members of the hierarchy of hyperviscous 
and G-hyperviscous equations. We have made a few tests analogous to those 
made in the case of the NS and GNS equations, and some of the results are 
briefly presented below. 

• (4) Local fluctuations. The theory of the fluctuations of a considered in 
item (2) above concerns the global phase space contraction rate, which in 
our models is extensive. Extensivity makes it very difficult to measure the 
fluctuations of a in systems with a large number of modes {i.e. many degrees 
of freedom). To overcome this difficulty, we define and study the "local phase 
space contraction" related to a portion of the volume occupied by the fiuid, 
where fluctuations are strongly enhanced. 

The investigation is important because in comparing fluctuations of phase 
space contraction between corresponding NS and GNS evolutions one must 
look at smaller subsystems. For instance, in the (truncated) NS equation 
there is no fluctuation of the total phase space contraction rate, because 
the friction coefficient is constant: the phase space contraction rate {i.e. the 
divergence of the equations of motion) is simply —2J2\]i\<N^^ if ^ is the 
cut off parameter. 

Furthermore, in real experiments, one is forced to look at the ffuctuations 
within small portions of fluid in order to obtain good statistics [15,16]. Our 
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results on this matter are only preliminary, but encouraging so that we 
report them in the Appendix. 

We integrated Eq.(5) by means of a pseudospectral code truncating it to a 
small number K of Fourier modes, Wk, of the velocity field. The computation 
of the entire Lyapunov spectrum is, generally, expensive (in terms of CPU 
time), and testing the validity of the FR (for the GNS equation) requires high 
statistics on the fluctuations of the phase space contraction rate a. Thus we 
have been forced to adopt a "low resolution" approximation of the dynamics in 
Fourier space: i.e. to take k, p, q G [— (A^ — 1), (A^ — 1)]^, where A^ is a small 
truncation value. In this way, the number of independent real components of 
the velocity field is 2K with 2K = 4N{N — 1) not counting the components 
corresponding to k = which, being constants of motion, have been fixed to 
be zero. We say that the equations are "truncated at K modes" calling "mode 
k" a (complex) Fourier component Wk of the velocity field. In our cases 2K 
ranges between 24 and 168: this corresponds to A^ — 1 between 2 and 6. Given 
our Reynolds numbers, K is much less than the value it should have (which we 
estimate to be about one order of magnitude higher) in order to interpret the 
results as relevant for developed turbulence (albeit in two dimensions), con- 
sistently with the theories of homogeneous turbulence. However, our present 
computer facilities make it impossible for us to reach much higher numbers of 
modes with the precision demanded by our analysis. 



3 Extension of equivalence between the NS and GNS dynamics to 
Lyapunov spectra 

The equivalence principle introduced in Refs. [11,21] yields quantitative rela- 
tions between certain observables, which can be tested in numerical exper- 
iments. The principle states that NS and GNS equations should have same 
statistics for the "local observables" 0, when the limit of large Reynolds num- 
ber is taken and the average phase space contraction of the GNS evolution 
equals the (constant) phase space contraction rate of the NS evolution. 

The phase space contraction rate associated with a differential equation is the 
divergence of the equations of motion: it can be straightforwardly computed 
from Eqs. (5), (6) to be a = (—1 + 2 J2\ki\<N 1) ■ a in the case £ = and m = 



i.e. for the time averages of functions of the velocity field which depend only on 
a finite number of its Fourier components or, more generally and more physically, 
which depend only on Fourier components of modes with k in the inertial range. 
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(which we call GED case), and 

\ i QiJ Qi 

in the GNS case, i.e. for £ = 1, m = (see Eq.(7))0. Denoting by 

a^^ = 2 \ki? = M (11) 

\k,\<N 

the (constant) phase space contraction rate of the NS evolution, and by 

= M-{a) + o{M) (12) 

the average phase space contraction rate of the GNS evolution, the equivalence 
conjecture studied in [22] can be stated as follows. 

Equivalence Conjecture (EC). The stationary probability distributions on 
the phase space associated with the NS equations and the GNS equations are 
equivalent in the limit of large Reynolds number, provided the value of Qq is 
chosen so that a^^ and (a^^^) coincide or, equivalently, provided the value 
of Qo is chosen so that (a) equals l|3 

The idea behind this conjecture is that if the fluctuations of a occur on time 
scales which are short compared with the macroscopic observation time scales, 
and if {a) = 1, the macroscopic observables should be equally well computed 
from the NS and the GNS dynamics. The large R limit is, then, required to 
ensure that the fluctuations of a do really occur on sufficiently short time 
scales. The results of [22] suggests that this is the case. 

To obtain such results, the steady state of the NS evolution was first deter- 
mined by following during a long time T the approach of the time average of Qq 
to its asymptotic value {Qo). At time T, the evolution was continued with the 
GNS dynamics with initial condition taken to be the vorticity field produced 
by the NS dynamics at time T. It was then observed that (a^^^) /cr^^ ^ (a), 
and that (a) differ little from 1 already at small (of order O(IO^) and even 
less) as long as the initial (fixed) value Qq^^ of Qq was close to its NS time 

^ In experiments on real fluids this quantity must be associated with some physical 

observable (a very delicate matter, [15,16]) 

^ i.e. equals the viscosity in our dimensionless units. 
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average (Qo)- Furthermore the details of the initial condition for the GNS 
evolution did not affect too much the statistical results produced by the GNS 
dynamics. In other words, it appeared that only the parameter Qq^^ needs 
to be adjusted to obtain equivalent NS and GNS evolutions. Although the 
meaning of "too much" and " close" was not sharply quantified in [22], the 
results suggested the validity of an EC also in the following alternative sense: 

Equivalence Conjecture (EC) . The stationary phase space probability dis- 
tributions of the NS equations and of the GNS equations are equivalent in the 
limit of large Reynolds number, provided the energy Qq is fixed in the GNS 
evolution to coincide with the average value (Qq) of the NS evolution at the 
same (large) R. 

In this paper, we first investigate the possibihty of extending the EC, EC' to 

the equivalence of the Lyapunov spectra. We evaluate the Lyapunov exponents 
of the NS equation by performing a long numerical simulation from time up 
to time T, yielding the NS Lyapunov spectrum and the average value (Qq), 
with an error bar SQq. From time T, we continue the NS evolution until Qo, 
which fluctuates, equals (Qo)- This usually takes a short time tg. Then, at 
time T + ts, we switch to the GNS evolution and compute the corresponding 
GNS Lyapunov spectrum. The reason to adopt this procedure, rather than 
just taking a random initial field with Qq^^ = (Qo) ^oy the GNS evolution, 
is that a priori the attractor might be not unique, and different initial fields 
may be evolved by the GNS dynamics towards attractors different from the 
one found by following the NS evolution. In other words, the NS and the GNS 
statistical properties may not always result comparable because of possible 
hysteresis phenomena. Our procedure was originally designed to avoid this 
problem, although it turned out subsequently that our choice of parameters 
yields no hysteresis, as discussed in Ref.[25]. 

The Lyapunov spectra produced by the NS and GNS dynamics in this way 
appear equivalent, meaning that the Lyapunov spectra of the NS and GNS 
dynamics appear to become identical for large R and K. At finite R and 
K, instead, the sums of the Lyapunov exponents should be related by (cf. 
eqs.(ll,12)) 

f:xr'-j:xr + o{M), where ^A^^-M. (13) 

1=1 1=1 1=1 

The typical situation is represented in Fig. 1 where the 2K — 2 nontrivial 
exponets are drawn. Note that two exponents are a priori known to vanish in 
the GNS case: one corresponding to the direction of motion and one to the 
direction of the conserved quantity; in the NS case one is a priori known to 
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Table 1 

Equivalence of NS and GNS dynamics, i.e. with I = \ and m = 0, for different 
Reynolds numbers. The last column gives the relative difference of the computed 
sums of the NS and GNS Lyapunov exponents (cf. eq.(13)). 





SQo/{Qo)ns 


Aa 


AQi 


o{M)/M 


408 


0.002 


0.039 


0.001 


0.060 


800 


0.002 


0.038 


0.001 


0.067 


2222 


0.003 


0.003 


0.002 


0.042 



Table 2 

Same as in Table 1 but for the ED and GED case with i = and m = 0. 

vanish, corresponding to the direction of motion. However in all cases the NS 
evolution shows another exponent very close to zero in agreement with the 
equivalence conjecture and we have taken it to be zero in drawing the figures 
1 and 2 below, see also footnote 11 below. 

The Lyapunov exponents are obtained by the method of Ref. [26] over times T 
of the order of up to 10^ Lyapunov time units.Q Error bars do not appear, since 
they are not larger than the size of the symbols in the figures, as explained 
below in Section 5. 

Consistently with the EC, and up to numerical errors, we observed that the 
viscosities and the average enstrophies of the NS and GNS evolutions with 
2K = 24 result equal if the average energies are equal. The numerical errors 
we refer to can be attributed to our approximate knowledge of the average of 
Qo (which is affected by an error bar SQo)', and to the finite values of R,N, 
and T (which, in turn, can only be reached through a sequence of time steps of 
finite size h). The corresponding deviations from the conjectured equalities of 
the viscosities and of the average enstrophies are measured by the terms Aa = 
{{a) - 1) and AQi = \{Qi)ns - {Qi)gns\/ {Qi)ns in Tables 1 and 2, where 
the subscripts NS and GNS indicate the two relevant kinds of dynamics. 



We define the Lyapunov time unit as 1/ Xmax > where Xmax is the largest Lyapunov 
exponent. 
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Fig. 1. Lyapunov spectra for the NS runs with normal viscosity {£ = 1) at = 2222 
(left) and = 5000 (right), and corresponding GNS runs with constrained energy 
Qq. The 2K — 2 nontrivial exponents are drawn by associating each vahic of the 
abscissa k = 1,2, . . . , K — 1 with the k-th largest exponent and the k-th smallest 
exponent A'^ = X2K-k-i- The symbols "+" refer to the NS spectra, the "*" to the 
GNS spectra, and the "x" to the sums (Afe + A^)/2 for the NS case. There is no 
"pairing" of the exponents to a common average value, unlike the cases of isokinetic 
Gaussian systems [27]. Essentially identical results have been obtained by studying 
the ED and GED equations. 

Further computations of the Lyapunov spectra were performed for 2K — 
168, starting from random initial data, evolving with the NS evolution, and 
then switching to a GNS evolution, as above. In this case, however, the fixed 
quantity Q„i of the GNS evolution was the energy Qq in one instance, the 
enstrophy Qi in another instance, and the palinstrophy Q2 in the last instance. 
The motivation of this computation was to check the validity of the conjecture, 
proposed and partially tested in [22] , that "several" pairs of equations of the 
hierarchies of hyperviscous and G-hyperviscous equations, with £ < m + 1, are 
equivalent to each other in the limit of large Reynolds numbers. 

If the EG can be extended to such general cases and to cover the Lyapunov 
spectra, the spectra should overlap (for large R) . Gonsidering the fact that the 
error bars SQo, SQi, SQ2 on (Qq), (Qi), {Q2) in these cases were larger than in 
the cases with 2K = 24 (the largest being SQi ^ 0.09{Qi)), our results for 
£ — 1 and m = 0,1, 2, seem to confirm this property. This is illustrated by Fig. 2 
which shows the equivalence, within our numerical accuracy, of the computed 
Lyapunov spectra of the NS and GNS systems. A similar equivalence was 
observed in several other cases which we do not describe here for brevity. This 
reinforces the previous results of [22], and provides a heuristic motivation for 
the validity of the conjectures EG, EG'. 

In Fig.l and Fig. 2 two exponents for the GNS equation vanish exactly reflect- 
ing that Qq is a constant of motion and neither expands nor contracts in the 
direction of the flow. There is a third exponent close to which, however, 
seems to be numerically distinct from 0. In the NS cases only one exponent 
is certainly (corresponding to the flow direction) but there are two other 
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^2K-k-1 
-2000 

20 40 k 60 80 

Fig. 2. Lyapunov exponents for one NS run with normal viscosity (+) at N = 7, 
B? = 10^, and forcing on modes (4, —3), (3, —4). The 2K—2 nontrivial exponents are 
drawn as in Fig. 1. The middle line (x) is the graph of (Afc+A^)/2. The other symbols 
refer to the corresponding GNS runs with fixed energy (*) i.e. m = 0,i = 1, fixed 
"enstrophy" (□) i.e. m = 1,£ = 1, ot fixed "palinstrophy" (■) i.e. m = 2,£ = 1. The 
spectra are obtained from runs of different lengths with T G [125,250], in units of 
iZ-^max) ^max being the largest Lyapunov exponent. The overlap of the four spectra 
(albeit approximate, because of the uncertianties SQo,6Qi, 6Q2) reflects the validity 
of the extension of the EC, EC' to the whole spectrum and to different members of 
the hierarchy of equations. The error bars, see Sec. 5, can be identified with the size 
of the symbols. 

exponents which are numerically very close to 0, see also footnote^^ below. 



4 Time reversal symmetry, cixiom C and the fluctuation theorem 

A fluctuation theorem (FT), [28], for the fluctuations of the phase space con- 
traction rate a'~^'^^{u) in our GNS truncated systems can be derived under 
the assumption that the Chaotic Hypothesis (CH) is verified by the GNS dy- 
namics at the considered Reynolds number, [14]. The CH can be formulated 
as follows: 

Chaotic Hypothesis (CH). A chaotic many-particle system or fluid in a 

stationary state can be regarded, for the purpose of computing macroscopic 
properties, as a smooth dynamical system with an attracting set on which the 
system is a smooth transitive Anosov system. 

This means that the attracting set can be considered, for practical purposes, 
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to be a connected "smooth surface" , on which the dynamics is hyperbohc and 
in fact verifies the Anosov property, see [12] for a mathematical definition. In 
other words, one supposes that the attracting set is a smooth surface possibly 
of dimension lower than that of phase space, thinking that its fractality is not 
important (here the difference between attracting set and attractor matters: 
the latter being very often fractal, cf. footnote 2 above). If the CH is valid 
for the GNS systems, the FT of [14] follows provided the attractor is dense 
in phase space and provided the evolution is time reversible. The FT imposes 
strong restrictions, which we illustrate below, on the form of the fluctuations 
of the contraction of the surface elements on the attracting set (cf. eq.(9)). 

Let t I— > SfU be the time evolution for the initial velocity field u generated 
by the equations of motion on phase space, and call T the total simulation 
time, which must be adequately longer than the characteristic time 6 of the 
fluctuations of a'^'^^ . Assume that relaxation to a statistically stationary state 
takes place in a time t^ <^T with t^^ 9 where ^ is a characteristic time for 
the evolution, e.g. 9 = l/Xmax, and introduce the infinite time average of the 
phase space contraction rate a'^'^^: 

(a^^^) = \im^ [ a^^^(Siu)rft = / /i(rfu)a^^^(5tu). (14) 





If T is large, we can subdivide the time interval [to, T], into a number of subin- 
tervals of length (r + ta), with r 9 and td a decorrelation time, separating 
the consecutive evolution segments of length r. Consider the dimensionless 
quantities 

^^^{z) = , I a^'''{Stu)dt, (15) 



tO + («-l)(T + td)+id 



with i = 1, (T—to) / [r+td) and arrange them in a histogram which allows us 
to approximate the probability distribution vr^ of the values ofp = a'^^^' . The 
total time T is empirically chosen large enough so that larger simulation times 
and different initial fields produce differences in the probability distributions 
that lie within the error bounds over a range of p's which is sufficient to 
perform meaningful tests. 

Within the range of values of p in which the number n of events per histogram 
bin is statistically meaningfully] and for sufficiently large r, it was observed 

^ Here, we chose the range with n > 10 (n being the number of events in a bin), 
and we weighted the contribution of each bin to the fit as inversely proportional to 
the associated error bar. In this way, the cases with low statistical relevance (e.g. 
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(cf. [22]) that the quantity 

^iP' ^) = ^ ^Jgns^ [log i^^iP)) - log (^^(-P))] (16) 



could be fitted by a function hnear in p with slope c(r), obtaining: 

^ [log (7r-(p)) - log (7r-(-p))] = c(r)p , (17) 



within the error bars. If the conditions of the CH and of time reversibility are 
verified, one can check whether or not c(r) converges to a limit Cqo, as r is 
increased. // the attractor was really dense in phase space^ then the validity 
of the FT would imply that Coo exists and, precisely, that Coo = 1- In such 
limiting case we refer to Eq.(17) by calling it the "fiuctuation relation" (FR). 

However, the attractors of systems with strong dissipation often are not dense 
in phase space. This is the case of the NS equations, in which the attracting sets 
have a large open complement. In such cases further assumptions are needed 
before the CH can be used to derive interesting predictions. We shall consider 
the possibility that there is proportionality between phase space contraction of 
volume elements around the attracting set and contraction of surface elements 
on the attracting set itself (which has dimension lower than that of phase 
space). This somewhat surprising possibility can be related to (and predicted 
from) other assumptions which involve geometric and dynamic hypotheses 
introduced in general in [24] and applied to the GNS equations in [11]. 

To understand the hypotheses, first note that smoothness of the attracting 
set and reversibility of the dynamics do not suffice to deduce Eq.(17). Indeed, 
the image under time reversal of the attracting set A will be a repelling set 
I A which, if the attractor is not dense in phase space, is disjoint from it, cf. 
Fig.3. 

Nevertheless the following may happen and be stable under perturbations of 
the equations of motion. Consider the stable manifold of the points on the 
attracting set: part of this manifold will be on the attracting set (which we 
suppose to be a smooth surface) and another part will stick out, cf. (a) in 
Fig.3. Likewise the unstable manifold of the points on the repelling set will be 
partly on it and partly stick out, cf. (b) in Fig.3. 

It is possible that the parts of the manifolds that stick out cross each other 
transversally along a manifold that connects the attracting and repelling sets. 



with n ~ 10) have little influence on the computed c(t)'s and, particularly, on its 
values extrapolated to large r. 

^ Hence the whole phase space would be an attracting set. 
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(a) 



(b) 



(c) 



A 



lA 



Fig. 3. Illustration of axiom-C systems. The bottom squares represent the attracting 
set A, while the top squares represent the repelling set I A. The leftmost figure with 
vertical segments represents a point x G A with a piece of its stable manifold 
connetcing A to I A. Similarly, the central figure represents a point x' G I A with 
a piece of its unstable manifold. The rightmost figure represents the corresponding 
intersection of the stable manifold of x with the unstable manifold of Ix, which 
associates a point Ix in I A with a point x in A. 

determining on each of them a single intersection point, cf. the line in (c) in 
Fig. 3. In this way, a correspondence / between attracting set and repelling 
set is constructed [24]. By construction the map / commutes with the time 
evolution maps Sf. The composition /* of / with the time reversal map / of 
Eq.(9) leaves the attractor invariant and anticommutes with time evolution 
restricted to it, i.e. 



A system verifying the above geometric properties, introduced and discussed at 
some length in[24] and in [29,12], is said to verify Axiom C. For such systems, 
time reversal symmetry can be "spontaneously broken" as a control parameter 
{R in our case) varies, meaning that the attracting set A becomes smaller than 
the whole phase space, and distinct from the repelling set I A. However, upon 
breaking, it spawns a new (lower) symmetry I* which acts on the attracting 
set and anticommutes with time evolution. In a sense, in such systems time 
reversal cannot he broken as it is always dynamically regenerated. 

In Refs. [24,29] this scenario was presented as the simplest possible way in 
which time reversal symmetry could be spontaneously broken via the splitting 
of the phase space of a reversible system into an attracting set, a repelling 
set and a large region in between, as in Fig. 3. Under this assumption, the 
CH holds for the dynamical system obtained by restricting the dynamics to 
the attracting set. To such system the FT apphes, and the fluctuations of the 
contraction rate of the area element of the attracting set obey the FR with an 
asymptotic slope Cqq = 1. 

A further difficulty is that, usually, it is not possible to measure the contrac- 
tion rate of the surface of the attracting set because this set is practically 



St ■ I* = I* ■ S_t on the attracting set. 
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inaccessible dynamical object. However such a measurement is not necessary 
if one can relate the area contraction on the surface of the attracting set to 
the total phase space contraction rate which is easily accessible, at least in 
numerical experiments. 

In Figs. 1,2 the 2K — 2 nontrivial Lyapunov exponents are arranged in pairs, 
each consisting of the fc-th largest exponent and the k-th smallest exponent 
(see [21] for a heuristic theoretical argument motivating this arrangement of 
our data). We imagine that the pairs with one positive and one negative 
exponent pertain to the restriction of the evolution to the attracting set.P^ 
The other pairs, consisting of two negative exponents, would then pertain to 
the directions of the stable manifold that stick out of the attracting set. 

Adopting the above view, we see that in a time T and on the attracting set the 
phase space contracts by J2*{^k,T + -^2K-i-k,T) where A^ y denotes the "local 
contraction" exponent of the surface of the attracting set over a time interval 
T around the phase point occupied by the system at the beginning of the 
interval. Here, J2* denotes the summation over the values of k corresponding 
to pairs of exponents with one positive and one negative exponent. The total 
contraction of phase space is Y^i i^k,T + ^2K-i-k,T) ■ The two contractions are 
asymptotic to T J2* + ^2K-i^k) and to TJ2k=ii^k + >^2K-i-k), respectively, 
if Afc denotes the k-th nontrivial Lyapunov exponent. Hence, their ratio tends 
to 



J2*{^k + ^2K-l-k 
J2k=li^k + ^2K-1- 



k 



J2k=l{^k + J^2K-l-k) ^ 



If the ratio C{T) is constant over the typical time scales over which we ob- 
serve the fluctuations of the phase space contraction, there is proportionality 
between the phase space contraction and the area contraction on the surface, 
and the proportionality constant is precisely c (which implies a weaker prop- 
erty than the one envisaged in[21]: the latter in fact seems too strong and not 
verified in the present case). We have tested whether this is the case, by tak- 
ing as time scale 1/Xmax, the inverse of the largest Lyapunov exponent. The 
result is that, after a short transient, the quantity C{T) is indeed practically 



By the axiom C, i.e. by the consequent existence of a time reversal symmetry on 
the attracting set, we know that on such set half of the exponents is positive and 
half is negative, once the vanishing exponent associated with the flow direction has 
been discarded together with the other vanishing exponent associated with the con- 
served quantity (energy, or enstrophy or any of the global quadratic quantities Qrn^ 
depending on the model). Arranging the remaining 2K — 2 exponents in decreasing 
order: Ai > A2 > ... > )^2K-2, one can form the pairs (Afc, A'^) of Figs. 1,2. 
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constant (with fluctuation amplitudes of a few percent of its average value) 
while the phase space contraction fluctuates with much larger amplitudes. 
This is also reflected, for instance, in the variances of the histograms of C{T) 
with short T"s for the cases we checked, which are more than one order of 
magnitutde smaller than the variances of the corresponding histograms of the 
phase space contraction. Therefore the FR with slope 1 for the (not directly 
accessible) surface contraction on the attracting set will imply a FR for the 
total (measurable) phase space contraction with a slope c given by Eq.(18). 

If in Eq.(16) we use the histograms for vr"^ obtained as outlined above, the 
result is well fitted by the linear law (17). This is the case starting from r's 
comparable to l/Xmax, and remains valid throughout the range of r's that we 
can consider. In this range, the values of c(r) can be computed through a least 
square fit and, as suggested by the CH, the results can be interpolated by a 
function like 

Coo + A/t , (19) 



where c^o is the asymptotic slope and A is a constant: in fact, in Anosov 
systems the limit of log7r'^(p) is reached with errors of order as r — > oo 
[28,14]. In all cases we considered, such an interpolation fits well our data. 
Figures 4-5 illustrate that the above statements are vahd even in the cases 
with small R^. The reason for the smallness of the r's that we can consider 
can be understood by looking at Fig.4 where, by comparing the graph on the 
left panel (smaller r) with the one on the right panel (larger r), we see that 
the range of p's over which Eq.(16) can be assessed rapidly decreases with r. 
This is due to the decrease of the frequency of the negative fiuctuations of 
fjGNS,T^ which rapidly becomes practically zero as r grows. 

Nevertheless, the interpolation by Cqo + ^/t, where the Coo and A are obtained 
from a least square fit, looks convincing within our range of r's and we there- 
fore investigate the consequences of assuming that it remains valid beyond this 
range. In this case, one can compare the values of Cqo with those of c obtained 
from Eq.(18). The results are reported in tables 3 and 4. One realizes that, 
within the accuracy of our numerical results, discussed more in detail in the 
next section, there is agreement between the predictions of Ref.[21] for c and 
our values for Coo- Of course, this is only a consistency check of the axiom C 
hypothesis which, as mentioned above, can only be indirectly tested in our 
cases. 

We also remark that the axiom C hypothesis leads to a simple expression for 

the dimension of the attracting set, according to which the ratio between the 
dimension of the attracting surface (on which the attractor is dense) and the 
dimension of the phase space {2K — 1 in our cases because of the existence of 
the conserved quantity Qq) is the ratio Dp — {2K — 1 — 2n_)/ {2K — 1), where 
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Fig. 4. Computed values of F{p;t), Eq.(16), for a GNS system with fixed energy 
{i = 1, m = 0) and i?^ = 686. The straight hues interpolate the data, and the error 
bars represent the statistical errors. The left panel has r 1.2, while the right panel 
has r 2.5, in 1/Xmax units. In order to obtain the fit in Fig.5 we repeated similar 
measurements for various values of r. Data refer to a 24-rcal modes truncation 
{2K = 24). The data for p < can be inferred from the symmetry p <-> —p. The 
kurtosis of tt'^ is defined by {{p — l)^)/((p — 1)^)^, and equals 3 if tt"^ is Gaussian. 
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Fig. 5. The value c(r) for the GNS run with fixed energy {£ = 1, m = 0), = 686. 
The times r, in the horizontal axis are given in 1/Xmax units. Data refer to a 24-real 
modes truncation {2K = 24) and to values of r to which already correpsonds a linear 
graph of F{p;t). The datum with smallest r has not been used in constructing of 
the fit (but it matches it, nevertheless). 

n_ is the number of pairs with two negative Lyapunov exponents. Obviously, 
the fractal dimension of the attractor cannot be larger than the dimension 

of the attracting surface. Therefore, under our assumptions on the nature of 
the attracting set, if the Kaplan- Yorke dimension is a measure of the fractal 
dimension of the attractor, its value divided by 2K — 1 should be < Dp. It is 
therefore of some interest to check whether this holds in our case or not. 
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1250 


0.639 


0.621 lb 0.024 


0.845 


0.652 (n_ = 4) 


2222 


0.834 


0.717 ± 0.011 


0.914 


0.826 (n_ = 2) 


5000 


0.842 


0.839 ± 0.015 


0.939 


0.826 (n_ = 2) 



Table 3 

Properties of the GNS dynamics with fixed energy = 1, m = 0). The value of c 
is obtained from Eq.(18) using the measured Lyapunov exponents. The value c^o is 
obtained by extrapolation of the data as in Fig. 5. The column labeled by KY, gives 
the Kaplan- Yorke dimension computed from our measured Lyapunov exponents. 
The column labeled by Dp is the ratio {2K — 1 — 2n-)/{2K — 1) where n_ is the 
number of pairs with two negative Lyapunov exponents. Data refer to a 24-real 
modes truncation (2K — 1 = 23). Our estimate of n_, given in brackets in the 
column of Dp, affects the calculations of both c and Dp. 

In most cases cases {e.g. in the cases of Fig.l) there are three exponents 
very close to and they are affected by a large relative error (see comments 
in sect. 5 below). We interpret as the two of smallest absolute value (as we 
know that two exponents necessarily vanish) . Although the signs of the second 
and third smallest exponents are stable with the length of our simulations, 
some uncertainty remains, at times, on which of the two approximates a zero 
exponent, and which approximates a finite exponent. For this reason, the 
values n_ reported in brackets in Tables 3 and 4 are uncertain by 1, and this 
implies an uncertainty also on the calcualtions of c and Dp (cf. point 3 of 
Section 6 below). 

Table 3 shows that in our system the axiom C prediction and the KY value 
do not seem really compatible. The comparison of the dimensions and its 
consequences has been suggested by F. Bonetto. The incompatibility between 
the axiom C and the KY methods for determining the fractal dimension of an 
attractor seems to be quite general and to apply to rather different kinds of 
chaotic motions. 

The same questions addressed here have been investigated in the case of the 
GED equations, obtaining the results reported in Table 4. 

It should be noted that the a priori pairing assumption which, together with 
axiom C, allowed us to interpret and organize our data, does not have firm 
theoretical grounds other than the discussion in [11]. Therefore the above 
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0.911 


0.823 ± 0.009 


0.943 


0.913 (n_ = 1) 


556 


0.813 


0.812 ± 0.006 


0.957 


0.826 (n_ = 2) 


800 


0.910 


0.910 lb 0.007 


0.983 


0.913 (ra_ = 1) 



Table 4 

Properties of the GED dynamics with fixed energy = 0, m = 0). The meaning of 
the columns is the same as in table 3 above. Data refer to a 24-reaI modes truncation 
{2K-\ = 23). 

analysis is highly conjectural and finds its main justification in its capability 
to fit the data. 

5 Technical coments on the simulations. Accuracy. 

The calculations summarized by Figs.l, 4-5, and by Tables 1-4, all refer to 
cases with rather small R and compared to the cases studied in Ref.[22]. 
Small R is needed to obtain values of Coo sufficiently small to be distinguished 
clearly from 1, while small N is needed to compute with reasonable accuracy 
the Lyapunov spectra. The actual values R and N that we chose constitute a 
compromise, but the result is that the frequency of the negative fluctuations, 
needed for a test of the fluctuation formula, is quite low, and very rapidly 
it becomes negligible for growing r's, even in long simulations. Given our 
present computer facilities, our standard simulation was of order of 10^ time 
steps, with sizes DT of order between 10~^ and 10~^, corresponding to times 
T of order between 10^ and 10^ in IjXmax time units. 

Such T's suffice for calculations of the average of quantities hke the Qm's, 
and to build with some accuracy the histograms of the fluctuations of a'^^^''^ 
for relatively small r's. The accuracy is reflected by the error bars placed on 
the computed quantities, like in Figs. 4-5, which cover the range of r's beyond 
which the negative fluctuations are too rare to have any statistical signiflcance 
(cf. footnote 8). These bars represent the statistical errors, obtained from the 
error bar on the histogram for tt"^, i.e. ±-y/n on a bin, if n is the number of 
counts in that bin. 

The times T are also sufficient for an accurate evaluation of the largest Lya- 
punov exponents. Indeed, out of 24, the largest (in absolute value) 20 expo- 
nents have an uncertainty ranging from 25% down to 0.06% of the computed 
value, for the largest one. The uncertainty on the remaining 4 (smallest) ex- 
ponents at times exceeds 100% of the computed value. The uncertainty on 
one exponent has been evaluated as the amplitude of the largest fluctuation 
of the computed value of that exponent in the second half of a simulation of 
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length T. In all cases, the error bars do not exceed the symbols size, if drawn 
in figures such as Fig.l. Nevertheless, at times, it is not obvious which two, 
of the three smallest exponents, vanish and which does not. This makes n_ 
uncertain by one unit. In particular, in the GNS case with = 2222 and 
in the GED case with = 408, could have been taken to be 3 and 2, 
respectively, rather than 2 and 1 as reported in Tables 3 and 4. In fact, the 
third and the second smallest exponents are both quite small and close in 
absolute value, but the third smallest is positive while the second smallest is 
negative. Therefore, considering zero the third smallest exponent, rather than 
the second, makes grow by 1. 

The slopes c(r), from which the values Coo of Tables 3 and 4 have been ex- 
trapolated, were computed by fitting our data to Eq.(17) and by using our 
histograms of a^^^''^ for vr"^. The limiting slopes Coo and its error bar reported 
in Tables 3 and 4 have been computed by a least square fit the function 
Coo + A/t to the c(r)'s. 



6 Concluding remarks 

1. The results on the equivalence between NS and GNS dynamics obtained in 
Ref. [22] have been strengthened through the analysis of the Lyapunov spectra 
of NS and GNS systems, albeit at very low spectral resolution. Similar results 
were obtained for spectral resolutions with up to 190 modes. 

2. The results of Ref. [22] on the validity of the FR for the GNS dynamics 
have been strengthened. Unfortunately, the relatively high dissipation of the 
systems considered here did not allow us to study the fiuctuations of ctGa^s.t^ 
with T » l/A^^^,, as desired, due to lack of appreciably frequent negative 
fiuctuations at such r's. Nevertheless, our data are well fitted by Coo + A/t, 
in the whole ranges of r we considered, in accord with the validity of the CH. 

3. The Lyapunov spectra have been used to compute the values c of Eq.(18) 
reported in Tables 3 and 4. These values equal the corresponding quantities 
Coo in agreement with the axiom-C hypothesis, except for the GNS case with 
_r2 = 2222 and for the GED case with R^ = 408. The uncertainty on c looks 
more pronounced than the uncertainty on Coo, because of the uncertainty on 
the number of positive Lyapunov exponents arising from the fact that often 
one of them is too close to 0.[^ The effect of this uncertainty on c is of about 

Two exponents have to vanish: the one corresponding to the conserved quantity 
Qm, and the one corresponding to the direction of the flow. However there are other 
exponents very close to zero, and one faces, therefore, the difficult task of distin- 
guishing numerically the small (positive or negative) exponents from the vanishing 
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9%, since it affects by 1 the number of pairs appearing in the mmierator of 
Eq.(18), and the contributions of the different pairs do not fluctuate much 
around their average value. In particular, in the GNS case with i?^ = 2222 
and in the GED case with R'^ — 408, taking n_ = 3 and n_ = 2, respectively, 
leads to a close agreement between c and Coo also in these two cases. 

4. Our codes, considered as dynamical systems in their own, seem to verify 
the axiom C, within the time scales of our simulations. Whether this property 
can be inferred to hold for the dynamics of real fluids is not clear, because of 
the low spectral resolution we have used. 

5. Global fluctuations of macroscopic systems can hardly be observed, al- 
though the experiments of Rcfs.[15.16] attempt at studying precisely this 
kind of fluctuations. However, as argued in the Appendix, the (observable) 
local fluctuations seem to obey laws similar to those of the global fluctuations 
discussed here. 

6. The extrapolation hypothesis Eq.(19) leads to a slope Coo which has been 
fltted to a value < 1. If conflrmed by further experiments, this implies that the 
stochastic fluctuation theorems which lead to Cqo = 1, such as those recently 
proposed in [30], would require extra ideas to be adapted to the description 
of the fluctuations of our deterministic dynamical systems. 



Appendix: Remcirks on local fluctuations and on spatio-temporal 
chaos. 

Recently, local versions of the FT have been considered to overcome the prob- 
lem that the large global fluctuations which are the object of the original FT 
cannot be observed in real macroscopic systems. One possible approach has 
been outlined in Ref.[31], where a fluctuation relation similar to (17) has been 
derived for the fluctuations of local quantities, in an infinite chain of weakly 
interacting chaotic maps, which is a paradigmatic system exhibiting spatio- 
temporal chaos. In particular, let Vq be a finite region of the chain centered 
at the origin, Tq > be a time interval, and define ?]+ and p as 

I To-l ^ To/2 

where V — Vq x Tq, and tiyq^x) is a properly chosen quantity which in [31] is 
called the entropy production rate in Vq and is identified with the phase space 

ones. 
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contraction rate in Vq. The following relation was obtained for the steady state 
probability distribution of p: 

ny{p) = e^(p)\y\+o(\dV\) ^ ^.^^ C(j^)-a-P) ^ ^ ^ p. ^ ^21) 

PV+ 

where \dV\ is the size of the boundary of V, p* > 1 and is analytic in p. 
This is perfectly analogous to the original FT, except for the boundary term 
\dV\ whose influence should diminish with growing V. 

The question of the vahdity of relations analogous to Eqs.(20,21) in 2-dimen- 
sional fluids is beyond the scope of the present paper. Nevertheless, it seems 
worthwhile to report our encouraging preliminary results on the validity of 
a local fluctuation relation for the NS dynamics. In fact, these results make 
plausible a linear relation with slopes smaller than 1, and indicate that in 
the NS case, in which dissipation occurs uniformly throughout the system, 
the local fluctuations may obey laws very similar to those obeyed by the 
global fluctuations considered in the previous sections. The interest lies in 
the possibility that, in agreement with the EC, the local fluctuations can be 
observed also in the NS dynamics, despite their irreversibility. 

The first study of a local fluctuation relation to obtain the value of Coo is in 
[32], and the problem has raised many discussions since. 

Let us consider the NS equation on a periodic cell of side 27r, and let Vq — 
[—L/2, L/2] X [— L/2,L/2] be a box with side L < 27r. We define the power 
dissipated in Vq by the forcing field f at time t as: 

PL{t) = J f{x) ■ u(x, t)Xy^ (x) c^x , (22) 

where Xv^-, characteristic function of Vq which, in spectral form, is: 

X^^(x) = Z)kXk(-^) e*^''. Defining local vorticity, local palinstrophy and a^, 
respectively, as 

Q,,,{t) = / (Va;)^(x, t) Xy, (x) dx , a,{t) = , (23) 

with a; = V X u, we propose to take rjvoi^) as proportional to the work per 
unit time performed by the external force within the region Vq divided by 
the corresponding local vorticity, i.e. we propose to take tiv„{'^) proportional 
to This is suggested by the facts that such a quantity has a physical 
meaning, and that the phase space contraction rate can in many systems 
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be identified with the entropy production rate. In fact, Qi.l can be called a 
"local temperature" being a quantity whose value measures the variations of 
the velocity field around its average. A possible alternative definition is a local 
version of Eq.(lO), hke (for instance) 

^^ = '^[[^) E - ^ + ' ' ■ (24) 

V^^^ \ki\<N Q^,lJ Qi,L 

However, several other local versions of Eq.(lO) could be defined: e.g. the last 
integral could be replaced by jV'^fu Xv^'^'^ by / / V^o; Xvb^'^- 

Since we approximate the various quantities introduced above by means of 
truncated Fourier expansions, L can only take the values ln/{N — 1) with 
/ = l,...,2(7V-l). 

Studying the behaviour of (Tl = ^l/ (o'l); being the infinite time average 
of (Tl, one realizes that a^, is essentially proportional to a^,, with proportional- 
ity constant equal to the ratio {ai)/ (ai), because the sum in (24) dominates 
over the other terms already at moderately large L. In particular, we per- 
formed a number of simulations with N = 7 (i.e. 84 complex Fourier modes) 
and = 2048, and we found that ai/io^L) and a^/ (o"^) are practically indis- 
tinguishable from each other for L > 57r/6. The statistics of the fiuctuations 
of (Jl was built from a run of 1.36 • 10^ time steps, corresponding to 2.4 • 10^ 
in 1/Xmax time units. 

As in the case of global fiuctuations (cf. Fig.4, and Ref. [22]), we found that the 
PDFs of the fiuctuations of the average of over times Tq are not Gaussian 
(cf. Table 5, and Fig. 7), and furthermore that the quantity 

F{p;V) = -—-Aog^^^ with P=Tf^^ I aLis)ds (25) 

is not linear in p at small V (i.e. small L, or small Tq), but it appears to be 
better and better approximated by a straight line of slope c{V) as V grows. 
Moreover, the slopes of the linear fits of the data appear to converge to values 
Cl for To — > cxo. Figure 6 illustrates the latter convergence at fixed box size 
L and growing Tq, showing that a linear fit of the data becomes acceptable 
at sufficiently large Tq, and suggesting that the slope of the fitting line may 
converge to a limit value. 

Figure 7 shows that the same holds if Tq is fixed and L grows, indicating that 
the nonlinear tails of F, may possibly be ascribed to non negligible boundary 
contributions, which become less important at larger L. In fact, we notice 
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Table 5 

Kiirtosis of Try for the cases of Fig.6. These values show that such PDFs are not 
Gaussian. 



that, for any V , F{p; V) is linear in a given neighborhood Iq oi p — 0, even 
if strong nonlinearities characterize F outside Iq. As V grows, it seems that 
the linear region persists, although it changes slope, while the nonlinear tails 
of F gradually disappear. Therefore, our data for F justify a test of the local 
fluctuation formula (21), in analogy with the theory of Ref.[31]. The result 
is summarized in Fig.8, where the slopes c of the linear fits of our data, and 
their variation with growing r are portrayed. There, we observe that our results 
for the local fluctuations are consistent with the EC and EC, and with linear 
relations with limit slopes smaller than 1. Slopes < 1 have also been previously 
found for local fluctuations in rather different systems, [32]. 

We conclude by remarking that: 

a. Our definition of the local phase space contraction rate (24) equals (up to 
a scaling factor) the physically measurable quantity «l for sufficiently large 
L, which in our case amounts to L = Ztt/G for 1 — 5,6, 12. 

b. The fluctuations of the local quantities pL and Pl/Qo,l, with Qq^l = 
J u^(x, t)x^^ (x) dx, behave in a perfectly analogous way as the fluctuations 
of aL. The latter remark is important because in the experimental work of 
Ref.[16] the quantity which is measured to attempt a check of the FR is es- 
sentially the work pL done by the external forces on a small region Vq per unit 
time, rather than the same quantity divided by an effective temperature. 

c. The fluctuation laws tested in this Appendix concern the {irreversible and 
constant friction) NS equation, and not the {reversible and fluctuating friction) 
GNS dynamics. This is consistent with the validity of the EC, EC. 



25 



Fig. 6. Numerically computed values of F{p;V), and corresponding error bars, for 
L = TT (left) and L = 5tt/3 (right). In each panel we have 5 curves, corresponding to 
To « 0.58,1.75,5.27,7.81,10.74 (in 1/A 

max units), ordered clockwise with growing 
Tq. Only three curves are clearly distinguishable because the three of highest Tq 
practically overlap in their range of coexistence, which shrinks with increasing Tq. 




1 p 2 1 P 2 



Fig. 7. Values of F{p;V), and corresponding error bars, for Tq ^ 3.7 (left) and 
To « 7.0 (right) in 1/A 

max units . In each panel we have 3 curves, corresponding to 
L = vr, 57r/3, llvr/G, and the values of k represent the kurtosis of the corresponding 
PDFs. 
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L___ . 
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Fig. 8. Slopes c{V) for fixed Vq as a function of Tq (in 1/Xmax time units). The left 
panel has L = 3tt/2 and the data are interpolated by c = 0.0025 + 0.024/To. The 
right panel has ll7r/6 and the data are interpolated by c = 0.0043 + 0.022/To. 
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